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Abstract 

The accumulation of small particles is analyzed in stationary flows through channels 
of variable width at small Reynolds number. The combined influence of pressure, 
viscous drag and thermal fluctuations is described by means of a Fokker-Planck 
equation for the particle density. It is shown that for extended spherical particles 
the shape of the fluid domain gives rise to inhomogeneous particle densities, thereby 
leading to particle accumulation and corresponding depletion. For extended spher- 
ical particles, conditions are specified that lead to inhomogeneous densities and 
consequently to regions with particle accumulation and depletion. 
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1 Introduction 



During recent years the investigation of microfmidic systems has increasingly 
attracted attention, being boosted by its promising and powerful so-called 
"lab-on-a-chip" [Tf2ll3] applications. A standard task that such a device should 
be able to perform is the separation of small objects immersed in a fluid 
according to specific properties of these objects like size or form. For this 
purpose, several mechanisms have been suggested. The proposed methods are 
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either based on downsized conventional laboratory tools like tubes, pumps and 
valves, or they employ effects that only function in the microfluidic regime. In 
combination with particular geometric constraints [HIqToT7]I8] and various ways 
of external forcing such as by electric fields [9~lll0|lll|ll2j . thermal gradients 



T3][Tj] or acoustic streaming induced by strong ultrasound waves 



thermal fluctuations play a key-role in many of these separation techniques 



Of great practical and principle interest are sorting methods that exclusively 
utilize inhomogeneities of flows in constrained geometries in combination with 
the omnipresent thermal noise. In particular, there are no external forces act- 
ing on the particle with these methods. The drift ratchet provides a typical 
example of this kind of device [T9f 201121 (22] . It consists of a large number of 
identical channels in a membrane, connecting two reservoirs both filled with 
water. The radius of each channel varies periodically in an asymmetric way. 
A carrier fluid, in most cases water, is periodically pumped to and fro. Al- 
though no water is transported on average, particles immersed in the water in 
general move to one side with a non- vanishing average velocity. Separation be- 
comes possible because the particle velocity and even its direction depend on 
the size of the particle [T9"ir2T)] . As in conventional ratchets [23.24:25(26 ^27112"%] 
the combined action of the periodic but asymmetric pore shape, breaking the 
left-right symmetry, the thermal fluctuations and the periodic pumping which 
drives the system out of thermal equilibrium are necessary ingredients for the 
observed effect [19|. 

The type of device studied in the present paper differs in several aspects from 
a drift ratchet. It does not require an oscillating flow, which is difficult to gen- 
erate and maintain experimentally [20|21|22j . but operates with a stationary 
flow field. The flow is confined to a channel with variable width - as for a drift 
ratchet - but with ends connected to each other in some way, leading e.g. to a 
ring-like or "figure eight" geometry, cf . Fig [TJ Spherical particles will then ac- 
cumulate in particular regions of the channel. In a previous experimental work 
the figure eight geometry was realized to separate particles immersed in a fluid 
by means of electric fields [29130] . In this experiment, a steady flow of the fluid 
was maintained by surface acoustic waves leading to acoustic streaming [15]. 

For the sake of simplicity we restrict our analysis to two spatial dimensions. 
The relevant mechanisms put forward in this work leading to an accumulation 
do not depend on the two-dimensional character and will work as well in 
three dimensions. In the experimental realization [2"9"f3"0] the flow channels are 
bounded by a rigid substrate only at the bottom, the remaining boundary is 
a free interface between water and air. A reduction of this very situation to 
two dimensions is not obvious. To still include the effects of different types 
of boundaries in the two-dimensional model we consider two cases, one with 
two no-slip boundaries corresponding to two impenetrable walls confining the 
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Fig. 1. (Color online) Sketch of the channel geometry. The middle part contains an 
asymmetric channel, which is doubly closed by two loops. The dashed line indicates 
the streamline connecting the stagnation points of the flow and separating the two 
flow chambers. The geometry may represent either closed channels or a wetting 
pattern on a flat substrate. The flow driving is exemplified to be acoustically driven 
as indicated by the two blocks of parallel lines representing surface-acoustic waves 
which propagate in the directions of the arrows [29|30| . 

fluid and the other one with a no-slip and a perfect slip boundary. The latter 
boundary mimics a free boundary with dominant surface tension such that 
the flow only leads to a negligible deformation of the equilibrium shape. 

The paper is organized as follows. In Section [2] we introduce the flow field 
which is then used for transporting particles. In Section [3] we analyze the 
hydrodynamic and random forces on the particles and formulate the long- 
time particle accumulation pattern in terms of a stationary Fokker-Planck 
equation. Section H] makes a distinction between possible accumulations caused 
by a volume and a boundary effect. In pressure-driven channel flows only the 
latter may occur, which is then investigated numerically in Section [51 



2 The flow fields 



We consider stationary incompressible flows in the limit of vanishing Reynolds 
numbers. The velocity v(x) and pressure p(x) fields are solutions of the sta- 
tionary Stokes equations, 

= -Vp(x) + 7?Av(x) + f (x), (1) 
= divv(x), (2) 

where r\ denotes the viscosity of the fluid and f is a conservative externally 
applied body force which gives rise to the flow fields p(x) and v(x). The 
velocity field is subject to the kinematic boundary condition, i. e. for immobile 
boundaries, 

v-N = 0, (3) 
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representing the fact that all boundaries of the fluid domain dQ are material 
lines of the flow. N denotes the normal vector. At immobile sticky walls, also 
the velocity components parallel to the tangent vectors vanish, 



v- 



T (a) = for all a. 



(4) 



For a free boundary, such as the interface between water and air, the tangential 
velocity at the boundary is determined by the mechanical stress balance. For 
the velocity field, this reduces to the perfect slip condition 



In general, at a free boundary the shape of the interface between two fluids 
depends itself on, and reacts back to, the flow and pressure fields in the fluid. 
Here we restrict ourselves to flows in a given domain on the boundaries of 
which the kinematic boundary condition ([3]) and either the no-slip @ or the 
perfect slip condition §5§ hold. The full description of free fluid surfaces with 
surface tension can be found in references [2H22]- In the numerical examples 
in Sec. |5]below, we will find that the type of the boundary condition has strong 
influence on the particle accumulation taking place near the boundary. 



3 Hydrodynamic particle transport 

The flow fields solving the Stokes equations ([T]) and §Z§ in combination with 
the boundary conditions ([3]) and (jlj) or ([5]) are used to advect small spherical 
particles. If such a particle is a point-particle at position X(t), it cannot be 
distinguished from the fluid material at this point. It is therefore transported 
with the velocity of the fluid itself, X(i) = v(X(t)). 

The motion of an extended spherical particle with small but non-vanishing 
radius R about X(i), however, is qualitatively different from that of a point- 
particle. The force on such a particle is the integral of the fluidic stress over 
the particle's oriented surface, denoted by F v (i) and F p (t) for the viscous and 
the pressure contributions of the stress, respectively. An additional random 
force F^(t) takes thermal fluctuations into account. For a small particle of 
mass m, the inertial force is negligible [33], resulting in an overdamped motion, 
described by the Langevin equation 



The pressure contribution F p (t) is obtained by expanding the pressure field in 
a Taylor series around the particle center X and integrating over the particle 




for all a. 



(5) 



= mX{t) = F v (t) + F p (t) + F 6 (t). 



(6) 
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surface S#(X), 
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p ~ 



- j> pNdA 



—R s VpW-0(R 5 ). 



(7) 



The viscous contribution F v (t) cannot be readily obtained, since the particle 
alters the velocity field: It poses a spherical no-slip boundary to the surround- 
ing fluid. In the Stokes regime, in which the acceleration field of the fluid 
caused by the movement of the particle can be neglected, we may employ 
Faxen's theorem of translational motion to describe this effect [Mf32] . The 
force F v (t) is then given in terms of a modified velocity field v, which is 
evaluated at the center of the particle, 



In contrast to the true velocity field v which describes the flow in the presence 
of the particle, the auxiliary field v is also defined within the region which 
is covered by the particle. Still, it contains all perturbations caused by the 
particle, cf. Appendix [S] for the definition and Ref. [32] for a detailed discus- 
sion. Unfortunately, determining v is as elaborate as solving the full stationary 
problem in the presence of the particle. In the next section we can, however, 
make use of the fact that v is a solenoidal vector field. 

The third contribution to the force on a particle in Eq. ([6]), i.e. the fluctuating 
force F^, reads at thermodynamic equilibrium [35 13 61137] . 



with the noise strength d = QirqRkvT ', and being Gaussian white noise, 



We neglect deviations from the equilibrium expression of the fluctuation- 
dissipation theorem, which become relevant for large velocity gradients only, 
see Refs. [38|39] , and references therein. Further, we assume a dilute solution 
of particles such that hydrodynamic particle-particle interactions can be ig- 
nored in the Langevin equation (jH]). The overdamped Langevin equation 
solved for the particle velocity is then equivalent to a Fokker-Planck equation 
for the probability density p(X, t) of the particle centers, 



F v = 6tt7]R\-X + v(X) + Av(X) R 2 /6 . 
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(9) 




(10) 
(11) 



dtp(x, t) = - div(p(x) d(x)) + J DAp(x). 



(12) 
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The diffusion constant is given by the Sutherland-Einstein relation 



d ksT 



(13) 



(67ir]R) 2 67n]R' 
and the drift velocity d reads, combining ([Tj) and (jSJ), 



R 2 

d(x) = v(x) + — Av(x) 



2R 2 



9i] 



V P (x) + 0(R 5 ). 



(14) 



The remaining term 0(R 5 ) comes from the Taylor expansion in Eq. ([7j). The 
viscous terms containing the modified velocity field v are exact. In the present 
paper, we focus on the long-time limit of the particle distribution, which is 
governed by the condition 9 t p(x, t) = 0, leading to the stationary Fokker- 
Planck equation 



4 Volume- versus boundary-accumulation 

As an accumulation of particles we refer to an inhomogeneous distribution 
of particles p(x). First, we note that the uniform distribution p(X) = const 
is a possible solution of the stationary Fokker-Planck equation (TT5]) in case 
of a solenoidal drift field, div(d) = [T2][T8] . Whether this solution is the 
physically realized one, depends only on the boundary conditions. Deviations 
from a uniform stationary particle density are caused either by particular 
boundary conditions which impose a non-vanishing gradient of the particle 
density at the boundary, or a non- vanishing divergence of the vector field d(x). 
We refer to the former and the latter case as boundary effect and volume effect, 
respectively. 

4-1 Volume effect 

The volume accumulation mechanism takes place for a drift velocity hav- 
ing a non- vanishing divergence, div(d) ^ 0. The homogeneous distribution 
p(X) = const then is not a solution of equation ffT5]) . irrespectively of any 
boundary condition. Whether the drift field is solenoidal, depends on the shape 
of the particle. A particle of complicated shape generally gives rise to a drift 
field with non-vanishing divergence. Examples can be found in the literature 



In the special case of a spherical particle, however, the drift velocity d(X) is 
given by equation (TT41 . Although the velocity field v(X) is not known explic- 
itly, it is known to be divergence- free. The same is true for its Laplacian. Only 




(15) 
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the pressure gradient may lead to a non- vanishing divergence of d(x). Taking 
the divergence of Eq. (fl4l) one obtains in leading order the Laplacian of the 
pressure p, which can be determined by the Stokes equation (pQ) to yield 



div(d) 



2R 2 



Ap 



2R 2 



div(f). 



(16) 



97/ 



9r/ 



In a long modulated channel, where the flow is driven either by a pressure 
difference, by a constant force, or by an inflow velocity profile, the drift ve- 
locity d(x) is indeed solenoidal. For these situations we state that all volume 
effects vanish exactly, and inhomogeneities in the particle density can only be 
caused by the boundaries. 

4-2 Boundary effect 

The boundary condition for the particle density in the fluid is given by the 
condition of vanishing normal particle flux, 



which guarantees that particles cannot cross the boundary. In case of a vanish- 
ing normal drift velocity N-d = only the uniform distribution p(x) = const is 
allowed. Vice versa, an inhomogeneous particle distribution can only emerge, 
if the normal drift at the boundary does not to vanish. 

The velocity field v, which is the leading-order term in the drift velocity (Tl4|) . 
satisfies the kinematic boundary condition ([3]) at the boundary of the fluid 
domain. As it does not provide a normal component, one might expect that 
neither a volume effect nor a boundary accumulation is achieved from the 
leading-order term in Eq. ( fl4l) . A close look on a spherical particle in the 
vicinity of a boundary, however, reveals that there is indeed such a normal 
component. Relevant in Faxen's theorem ([8]) is the modified velocity field v 
at the center of the particle, which cannot touch the physical boundary but 
stays at least one radius apart from it. The no-flux boundary condition (TlTl) 
must therefore be imposed on an effective boundary dVt' which has constant 
distance R from the physical boundary dfl of the fluid. Both are sketched in 
Fig. [2j Neither the streamlines of the flow nor the field-lines of the particle drift 
velocity need to have a constant distance from the boundary. Instead, some of 
the field-lines cross the boundary dQ' of the region which is accessible to the 
particle centers. In Fig. [2] one of these crossing points is marked. There, the 
normal component of d(X) is positive, i.e. pointing to the boundary whence 
leading to a deposition of particles at the boundary. Accordingly, a negative 
normal component results in a depletion zone. Both effects will be found in 
the numerical simulations below, see Fig. [5] 



= N-(pd-DVp), 



(17) 
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Fig. 2. (Color online) The boundaries of the fluid domain and of the accessible 
region f2' for the particle center. Their distance equals the particle radius R every- 
where. At the boundary <9f2' the normal component of the drift velocity, and thus 
also of the drag force, are not zero. This is exemplified at the point marked by a 
cross. 

The given argument in the foregoing paragraph is a geometrical explanation for 
the leading-order terms of the hydrodynamic interaction between particle and 
wall. A general theory of hydrodynamic interactions in asymmetric geometries 
is presently not available. A known effect of the hydrodynamic interaction 
is a decrease of the fluid velocity induced by the presence of the particle. 
Additionally, one obtains the Saffman lift force [15] which drives the particle 
away from the boundary. The magnitude of this force is proportional to the 
squared particle radius and to the mismatch between the true particle velocity 
and the unperturbed fluid velocity. Reference [16] provides explicit expressions 
for the correction terms from hydrodynamic interaction between particle and 
boundary. The force densities given there comprise a component, diverging at 
the effective boundary dQ', corresponding to the hard-core interaction. The 
next leading-order term in the particle radius R is the random force, being 
proportional to R 1 / 2 . The tangential force leading to the slowing down of the 
particle is proportional to R, and the Saffman lift force scale as R 2 . 

About the normal projection of the second-order terms in the drift veloc- 
ity d(x), being proportional to R 2 in Eq. (j!4p . we cannot say much. Gener- 
ally, they may give rise to a boundary effect. The magnitude of this effect, 
however, scales quadratically with the particle radius R and is therefore ex- 
pected to be smaller than the hard-sphere boundary effect. Furthermore, for 
pressure-driven flows in long channels, the pressure gradient at the boundary 
is oriented along the channel rather than normal to it. 



5 Numerical examples in modulated channels 

The aim of this paper is to provide a qualitative picture of the influence of 
the domain shape on the particle density. In order to achieve this goal we 
have to restrict the following numerical analysis to the leading-order terms of 
the forces. These are the geometrical interaction between particles and wall, 



as depicted in Fig. [2j and the random forces. By omitting terms which are 
quadratic in the particle radius, we may also replace the unknown velocity 
field v by the unperturbed velocity v. 

Using these simplifications, the full interaction between particle and wall is 
replaced by a hard-core interaction preventing the particle from penetrating 
the wall. We are aware that this approximation may be quite crude for the 
quantitative description of a real particle in a real fluid. However, the analysis 
above shows that the only terms which lead to an accumulation are the normal 
forces on the particles in the vicinity of the boundary. In the hard-core model, 
this happens instantaneously as the particle touches the wall. The force is 
localized on the effective boundary dfl'. Corrections to the force, as detailed 
in Ref. [16] smear out the force such that the particle is slowed down already in 
the vicinity of the wall. This will not alter the qualitative picture obtained here 
by the hard-core interaction. Whether this approximation is really justifiable, 
can be decided only by extensive numerical calculations of the hydrodynamic 
interactions, a task which is beyond the scope of the present study. 

Next, we solve the Fokker-Planck equation 

= -vVp + L>Ap in ft', (18) 

with the boundary condition 

= N'-(pv-DVp) on dVL'. (19) 

The vector N' denotes the normal vector of the effective boundary dVt' . The 
velocity field v is taken from the numerical solution of the Stokes equations 
in Q with the boundary conditions from Sec. [21 the domain Q being a periodic 
two-dimensional channel of modulated width. Its two boundaries have different 
geometries, one is straight, while the other is curved, see Fig. HI We choose the 
shape of the curved boundary to be given by the periodic function g, which is 
implicitly defined by the relation 

g(z) = sinfonz - a g{zj) . (20) 

The number a parametrizes the asymmetry of the shape, yielding a sinusoidal 
shape for a = 0, while shapes with a/0 have alternating steep and flat flanks. 
For a > the flat flanks have positive and the steep flanks have negative incli- 
nation, and vice versa for a < 0. The values g(z) can be evaluated iteratively, 
starting with sm(2iTz). 

The parameter regime of equation (fT8{) is characterized by the dimensionless 
Peclet number, which expresses the ratio of advective to diffusive transport. 



9 



Fig. 3. (Color online) A detail of the computational mesh near the boundary. The 
outermost layer of finite elements represents the zone which cannot be entered by 
the particle centers. The element sides of this layer constitute the two boundaries 
<90 and dVt' , compare with Fig. [2J 



Velocity magnitude f(x) in 10 




Fig. 4. (Color online) Example-flows in periodically continued two-dimensional chan- 
nels, modulated according to Eq. (|2U|) with asymmetry parameter a = 0.8. In 
panel (a) both boundary conditions are no-slip, while in (b) the straight bound- 
ary carries a perfect-slip condition. The flow direction is from left to right, driven 
by a pressure difference of Ap = 1 Pa per length L = 10~ 4 m of the unit cell. The 
viscosity is that of water. The black lines represent streamlines. 

We here define the global Peclet number as 



Pe 



xv 



3ttL 2 R Ap 
128 k B ~f~' 



:2V, 



For the stationary particle distribution, there is no difference between a weakly 
driven system and one at a high temperature, as long as the ratio Ap/T is 
the same. In the numerical calculations for several different temperatures, we 
therefore used the same numerically obtained velocity field, driven by a unit 
pressure difference Ap = 1 Pa along the channel length L = 10 _4 m. 



The computational mesh for the numerical solutions has to provide both 
boundaries, dQ and dQ'. Figure [3] shows a detail of the employed mesh near 
a boundary. The outermost layer of finite elements build the zone of constant 
width R which cannot be entered by the particle centers. This zone can be 
identified in Fig. [5] as the zone near the boundary, showing a vanishing particle 
density. 
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Particle density p(x) 

-1.022 

1.015 
-1.009 

1.003 
■0.996 
■0.990 

Particle density p(x) 

-1.080 
■1.050 
-1.030 
1.000 
■0.976 
■0.950 

Fig. 5. (Color online) Example-densities of particle centers in the flows of Fig. 01 
normalized to an average value of unity. The densities have been calculated as the 
solution of the Fokker-Planck equation (|15p with Ap/T = 3.3 x 10 -6 Pa/K which 
corresponds to a Peclet number of Pe ~ 16. At the boundaries, the black zone with 
zero density indicates the area which the particle centers cannot enter. It is one 
particle radius wide, here R = 10~ 2 L. The thin black lines indicate level lines of the 
density. 

5.1 Longitudinal accumulation 

Figure H] presents two velocity fields in the unit cell of a channel with asymme- 
try parameter a = 0.8. They differ by the boundary condition at the straight 
wall: Panel (a) presents no-slip and (b) perfect slip boundary conditions. The 
corresponding solutions of the stationary Fokker-Planck equation (|18|) with 
boundary condition (|T9|) are illustrated in Fig. [5j The qualitative form of the 
stationary distributions is quite intuitive. Left of the bottleneck we find an ac- 
cumulation zone of the particles. Resulting from the main flow direction from 
left to right, the particles are concentrated at the left side, just upstream in 
front of the bottleneck. This result can be understood by means of Fig. [2J Left 
and right of the bottleneck we find regions where streamlines cross the effective 
boundary dQ', leading to a non- vanishing drift component in normal direc- 
tion. This normal component results in the accumulation that we see in Fig. [5] 
for Pe ~ 16. This corresponds to a slightly advection-dominated transport of 
particles. For smaller Peclet numbers, the diffusion becomes more dominant, 
resulting in a smoother particle distribution than in Fig. [5j For larger Peclet 
numbers, the maxima of the distribution near the boundaries become more 
pronounced. 

5.2 Influence of the boundary conditions 

Comparing the particle densities in Fig. [5] for different boundary conditions of 
the velocity field, one finds the extremal values to be more pronounced in the 
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case with a perfect-slip condition. The reason for this difference is that the slip 
condition allows a larger velocity at the boundary, giving rise also to a larger 
normal component, which then causes the accumulation effect. We expect 
that in general free surfaces, providing perfect slip boundary conditions, cause 
larger accumulation effects than sticky walls, as long as the particles do not 
leave the fluid nor cause major deformations of the free surface. 



5.3 Perpendicular accumulation and sorting 



The two accumulation patterns in Fig. [5] do not only differ with respect to 
their magnitude. So far, we have only looked at the inhomogeneity along the 
channel orientation, by stating that there are concentration and depletion 
zones in front of the bottleneck and behind, respectively. For a small channel 
as in experimental applications also the distance of the accumulation centers 
will be prohibitively small. A performance characteristic of much more interest 
is the accumulation perpendicular to the channel orientation. This may be 
achieved by a branching channel where the flow is split along the streamline 
that ends in a stagnation point. Such a structure has been experimentally 
realized as a eight-shaped geometry and used for separating particles with 
external electric forces [29]. Particles immersed in this flow do not generally 
follow the streamlines and thus can cross the separating streamline and pass 
from one basin into the other. This leads to a relative accumulation of particles 
in one of the basins, if they were originally equally distributed. 



The quality of separation is characterized by the amounts of particles that 
are above and below the separating streamline indicated in Fig. |H For this 
purpose we separately integrate the particle density in the two regions Q' + and 
ft'_, representing the accessible regions above (+) and below (— ) the separating 
streamline, respectively. The resulting two probabilities per area P + / Q' + and 

P_/ ft'_ of finding a particle above or below the streamline can be used to 
provide a measure for the relative accumulation of particles in one of the 
basins. We use the quantity 



P_ 



r :- 



ft'. 



ft'_ 



(22) 



as a measure for the relative accumulation in the upper basin, which is the 



one with the curved boundary. The denominators 



ft' 



and 



ft' 



represent the 



volumes of the regions above and below the separating streamline, respectively. 



Figure [6] displays the resulting relative accumulation for flows with two no- 
slip boundaries in three different channel geometries, indicated by different 
values of the asymmetry parameter a. The upper panel shows the ratio r as 
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Fig. 6. (Color online) The relative accumulation r perpendicular to the main flow 
direction for flows with no-slip boundary conditions. Values smaller than unity cor- 
respond to an accumulation of particles below the separating streamline, near the 
straight wall. In panel (a) the temperature is varied for a fixed particle radius 
R = 10~ 2 L. In panels (b) and (c) the particle radius is varied for two selected 
temperatures, Ap/T = 3.3 x 10 -7 and 3.3 x 1CT 8 , respectively. 

a function of the inverse temperature. The first observation is that the effect 
vanishes for vanishing driving - or equivalently, for infinite temperature. For 
the majority of the parameter values the result is smaller than unity. This cor- 
responds to an accumulation of particles at the side of the straight wall. This 
appears as a general tendency, which was found also for other shapes. The 
relative accumulation effect is at most one percent for the smallest tempera- 
ture that was used in the calculation. As expected, it vanishes for very small 
driving strengths Ap/T. In order to demonstrate that this small effect is not 
an artifact of the numerical calculation, also the flows in the mirrored channels 
with inverted pressure differences were considered. For symmetry reasons both 
configurations yield identical accumulation ratios. The numerical differences 
between these are found to be smaller than the line thickness in Fig. For two 
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temperature values, also the particle radius was varied. Panels (b) and (c) of 
Fig. [6] depict the ratio r as a function of the radius. Again, the effect vanishes 
with vanishing radius. This is the expected behavior, because a point-particle 
can come arbitrarily close to the physical boundary. 

An interesting aspect of the accumulation results in Fig. O is the fact that 
the channel with a = 0.8 behaves differently from the one with a = —0.8. 
We find as a general tendency that the channel which suddenly widens and 
slowly narrows (a = —0.8) in flow direction yields better accumulations than 
the suddenly narrowing one. The channel with the sinusoidal shape (a = 0) 
yields accumulation results somewhere in between. 



5.4 Accumulation inversion 



Another remarkable property of the accumulation mechanism can be observed 
for the flow in the channel with asymmetry parameter a = 0.8. Here, we find 
values of r also being larger than unity. Particles in this parameter regime are 
transported towards the curved boundary rather than towards the straight 
one. The occurrence of both values corresponds to an inversion of the transport 
direction. Panel [6fc confirms that the values larger than unity, which we found 
in[HK persist also for several smaller radii. Note that for each radius, a different 
numerical mesh was used in order to represent the correct boundary dQ'. The 
inversion found here is not readily usable for sorting particles because it takes 
place upon varying the parameter Ap/T and not the radius R. Moreover, the 
effect seems far too small to be of experimental relevance. Still, the occurrence 
of accumulation inversion is an interesting effect. 

Employing a perfect-slip boundary condition at the straight wall leads to a 
qualitatively different parameter dependence of the accumulation. Fig. [7J de- 
picts the relative accumulation for the same driving parameters as Fig. [6j 
Again, the accumulation vanishes for vanishing driving and for vanishing par- 
ticle radius. The direction of the accumulation depends now much stronger on 
the driving strength Ap/T. As a qualitative picture, we here observe a par- 
ticle of a given size to be pushed towards the curved or towards the straight 
wall, depending on the pressure difference. This accumulation inversion can 
also be found in Fig. [7b as a function of the particle radius. Thus, we may 
find particles with different radii accumulated in the two different parts of the 
channel. For small particles this effect though is small. 
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Fig. 7. (Color online) The same as Fig. [6] for a flow with perfect-slip boundary 
condition at the straight wall. 

6 Conclusions 

We analyzed the accumulation of small spherical particles in two-dimensional 
fluid flows confined to channels with one curved boundary. It turned out that 
accumulation due to volume effects is impossible if the flow is driven by a 
pressure difference or a homogeneous force. In this case, only boundary ac- 
cumulation effects are possible. For that to happen the normal projection of 
the particle drift at the boundary must not vanish. Otherwise, only the trivial 
uniform distribution of particles is achieved. In order to simplify the numer- 
ical calculations we employed a hard-core interaction model for the intricate 
interaction between a particle and a wall. Both imply non-vanishing normal 
components of the drift velocity near the boundary, thus giving rise to a small 
accumulation effect, which scales with the particle radius. 

We compared no-slip and perfect-slip boundary conditions for the velocity 
fields in the channels. The latter is a simplification of the free-surface boundary 
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condition which is realized in some experiments. We observed that the type 
of the boundary condition on the flow has severe implications on the resulting 
accumulation pattern of immersed particles. Free surfaces generally caused 
larger accumulation effects showing also qualitatively different and interesting 
properties, such as an inversion of the accumulation direction, depending on 
the radius of the considered particle. This inversion might be put to beneficial 
use for a separation of particles of varying size. 
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A Velocity fields 



The following three velocity fields are relevant in Faxen's theorem to express 
the motion of a spherical particle in an external velocity field: First, in ab- 
sence of the particle, there is the unperturbed velocity field solving the Stokes 
equations ^ and (J2]). The flow field satisfies the integro-differential equation 



v k(y) = <j> dA(x) Nj(-x) K. ik (x - y) cr y (x) - T ijk (x - y) i>j(x) , (A.l) 
an 

where the external conservative force is taken into account by defining a mod- 
ified pressure containing the potential of the force. The kernels K, and T may 
be chosen to be the Green functions for the unbounded Stokes problem, 

87r?7 V r r 6 J 
= ~ rJ p- (A.3) 

The transformation of the Stokes equations into an equation on the boundary 
only is essential for Faxen's theorem. This implies that only conservative forces 
may drive the fluid. 

The second flow field which is important is the true velocity field in presence 
of the particle. We denote it with v. It satisfies a similar relation as (lA.lj) . 
but the particle cuts a spherical region out of the domain fl. The boundary 
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integral is then carried out also over the surface of the particle, 



S fc (y) = jdA(x)N 3 (x) 

anus R (x) 



/Q fc (x - y) try(x) - T ljk (x - y) ^(x) 



(A.4) 



The third relevant velocity field is the one that is used in Faxen's theorem. It 
must be evaluated at the center of the sphere. This effective velocity field v 
is obtained by carrying out the integral in ( 1A.4j) over the particle surface 
explicitly as 



(y) := j> <L4(x) Nj(x) [/C ife (x - y) a^(x) 
an 



- T ijk (x - y) 



vAjc) 



(A.5) 



For the details of this calculation see Ref. [32J. This effective velocity field is 
is then defined in the whole domain Q and can be evaluated at the center of 
the particle. It is used for expressing the force on a particle in Eq. (jSj). Note 
that the boundary values of the true velocity field v enter the definition of 
v. As a solution of the Stokes equation with appropriate pressure field, it is 
solenoidal. We make use of this property in Sec. 14.11 
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